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Abstract. In this paper we continue the study started in [IB]. We consider a planar, 
bounded, m-connected region J7, and let dfl be its boundary. Let T be a cellular decompo- 
sition of £7U<9f2, where each 2-cell is either a triangle or a quadrilateral. From these data and 
a conductance function we construct a canonical pair (S, /) where S is a special type of a 
(possibly immersed) genus (m— 1) singular flat surface, tiled by rectangles and / is an energy 
preserving mapping from onto S. In [16] the solution of a Dirichlct problem defined 
on T' - 1 was utilized, in this paper we employ the solution of a mixed Dirichlet-Neumann 
problem. 



0. Introduction 

Before stating our main result, we need to define a special kind of two dimensional ob- 
jects, surfaces with propellors. A flat, genus zero compact surface with m > 2 boundary 
components endowed with conical singularities, will be called a ladder of singular pairs of 
pants. A sliced Euclidean rectangle is a Euclidean rectangle in which two adjacent vertices 
are identified, and possibly a finite number of points on the opposite edge have been pinched. 

Definition 0.1. A singular flat (possibly immersed), genus zero compact surface with m > 2 
boundary components having conical singularities, will be called a surface with propellors, 
if it has a decomposition into the following pieces: a ladder of singular pairs of pants, sliced 
Euclidean rectangles, Euclidean rectangles and straight Euclidean cylinders. 



We consider (as in [16]) a planar, bounded, m-connected region f2, and let <9f2 be its 
boundary. Let dQ = E\ U E 2 , where E\ is the outermost component of dQ. Henceforth, we 
will let {«!, . . . , ai} be a collection of closed disjoint arcs contained in E\, and let . . . , f3 m } 
be a collection of closed disjoint arcs contained in E 2 . Let T be a cellular decomposition 
of Q U dQ, where each 2-cell is either a triangle or a quadrilateral. Invoke a conductance 
function on , making it a finite network, and use it to define a combinatorial Laplacian 
A on T(°\ These data will be called Dirichlet-Neumann data for (VL,dVL,T). Let k be a 
positive constant, and let g be the solution of a mixed Dirichlet-Neumann boundary value 
problem (DN-BVP) defined on and determined by requiring that 

(!) g\ ai = k, for alU = 1, . . . , /, and g\ E2 \{p 1 \j...up m ) = 0, 
da 

(2) Q^lft = 0, for all j = 1, ... , m, 

(3) Ag = at every interior vertex of T'°\ i.e. g is combinatorially harmonic, and 

(4) E. 680 fS(fln)(*) = o, 
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where (4) is a necessary consistent condition. Let E(g) denote the Dirichlet energy of g. We 
may now state the main result of this paper: 

Theorem 0.2. (Main result) Let (Q,dQ,T) be a bounded, m-connected, planar region en- 
dowed with a Dirichlet-Neumann data, with m > 2. Then there exists a surface with propel- 
lors Sq, and a mapping f which associates to each edge in a unique Euclidean rectangle 
in Sq in such a way that the collection of these rectangles forms a tiling of Sq. Furthermore, 
f is boundary preserving, and f is energy preserving in the sense that E(g) = Area(Sn). 

Throughout this paper a Euclidean rectangle will denote the image under an isometry 
of a planar Euclidean rectangle. For instance, some of the image rectangles that we will 
construct embed in a flat Euclidean cylinder. These will further be glued in a way that will 
not distort the Euclidean structure (see $2] and §4] for the details). 

In our setting, boundary preserving means that the rectangle associated to an edge [u, v] in 
T^ x ' with u G dfl has one of its edges on a corresponding boundary component of the surface. 
In the course of the proof of Theorem l0.2l it will become apparent that the number of singular 
points and their cone angles, the lengths of shortest geodesies between boundary curves in 
the ladder and the number of propellors, may be explicitly determined. In particular, the 
cone angles obtained by our construction are always tt/2, or even multiples of n. Some classes 
of such surfaces are called translation surfaces, and for excellent accounts see for instance 
HE], EI] and [23]. 

Also, the dimensions of each rectangle are determined by the given DN-BVP problem on 
T^°\ Concretely, for [u,v] G T^ l \ the associated rectangle will have its height equals to 
(g(u) — g(v)) and its width equals to c(u,v)(g(u) — g{v)), when g(u) > g{v). Some of the 
rectangles are not embedded. We will comment on this point (which is also transparent in 
the proof of the theorem above) in Remark 15.121 In a snapshot, some of the rectangles which 
arise from intersection of edges with singular level curves of the DN-BVP solution are not 
embedded. 

A surface with dents and pillows will denote the surface obtained by doubling a surface with 
propellors along its boundary. The following Corollary is straightforward, thus establishing 
the statement in the abstract of this paper. 

Corollary 0.3. Under the assumptions of Theorem \0.^\ there exists a canonical pair (S, f ), 
where S is a flat surface with dents and pillows of genus (m — 1), having conical singularities 
and tiled by Euclidean rectangles. The mapping f is energy preserving from into S, in 
the sense that 2E(g) = Area(S'). 

Proof. Given (f2, <9f2, T), glue together two copies of Sn (their existence is guaranteed by 

Theorem 10.21) along corresponding boundary components. This results in a flat surface 

S = Sn U Sq of genus (m — 1) and a mapping / which restricts to / on each copy. 
an 

□ 

The following two theorems are foundational and serve as building blocks in the proofs of 
the theorems above. In [161 Theorem 0.4] we proved the following: 

Theorem 0.4. (Discrete uniformization of an annulus [8]) Let A be an annulus and let 
(Q,dQ,T) = (A,dA — E\ U i?2,T). Let Sa be a straight Euclidean cylinder with height 
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H = k and circumference 

C=£g(x). (0.5) 

Then there exists a mapping f which associates to each edge in a unique embedded 
Euclidean rectangle in £4 in such a way that the collection of these rectangles forms a tiling 
of Sa- Furthermore, f is boundary preserving, and f is energy preserving in the sense that 
E(g)=Avea(S A ). 

A topological planar closed disk with four distinguished points on its boundary, its corners, 
will be called a quadrilateral. Let 1Z be a quadrilateral endowed with a cellular decompo- 
sition. Let &1Z = d7Zb ottom U dTZ\ c f t U dlZ top U 97^ ri ght be a decomposition of dlZ into four 
non-trivial arcs of the cellular decomposition with disjoint interiors, in cyclic order. If the 
intersection of any two of these arcs is not empty, then it consists of a corner (all of which 
are vertices). A corner belongs to one and only one of the arcs. We solve a mixed Dirichlet- 
Neumann boundary value problem with E\ = E 2 , ot\ = dTZ T i g ht, Pi = <97^bottom , P2 = d7Z top 
(and hence that and g\dTz Mt — 0). The second foundational theorem is proved in this paper: 

Theorem 0.6. (Discrete uniformization of a rectangle [8]) Let 1Z be a quadrilateral. Let Sr. 
be a Euclidean rectangle with width W = k and height 

H = £ §^)(*)- (°- ? ) 

Then there exists a mapping f which associates to each edge in a unique embedded 
Euclidean rectangle in S-r in such a way that the collection of these rectangles forms a tiling 
of Sr. Furthermore, f is boundary preserving, and f is energy preserving in the sense that 
£(<?)= Area(^). 

Given (Q,dQ,T), we will work with the natural affine structure induced by the cellu- 
lar decomposition. Let us denote this cell complex endowed with its affine structure by 
C(tt,dQ, T). Next, study the level curves of g on a 2-dimensional complex which is ho- 
motopically equivalent to C(f2, <9f2, T), embedded in R 3 , obtained by using g as a height 
function on (Q,dtt, T). We will work with the level curves of g or equivalently, with their 
projection on C(Q, dfl, T). 

Once Theorem 10.61 is proved, we proceed to prove Theorem IU.2I as follows. We first con- 
struct a topological decomposition of Q. For i = 0, . . . , k, consider the sub-domain of Q 
defined by 

fli = {x e n I Pi < g(x) < p i+1 }, (0.8) 
(where the value at x which is not a vertex is defined by the affine extension of g). In 
general fij is multi-connected and (by definition) contains no singular vertices in its interior. 
Let g.i = glfy be the restriction of g to U <9f2j. The definition of gi involves (as in the 
proof of Theorem 10.61) new vertices (type I), and new edges and their conductance constants. 
In particular, each ^ is the solution of a D-BVP (Dirichlet boundary value problem) or a 
DN-BVP, on each one of the components of fij. 

By applying a topogical-combinatorial index lemma (Lemma 13. 7j) and a splitting argum- 
net, we will conclude that Q may be decomposed into a union (with disjoint interiors) of 
annuli, quadrilaterals or sliced quadrilaterals. We will finish the proof by showing that the 
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gluing is geometric, i.e. that for all i, any part of dfli has the same flux-gradient metric 
(Definition ll.6p . with respect to the boundary value problems induced on each of the two 
components it belongs to. 

Throughout this paper we will assume that the reader is familiar with the results and 
terminology of [16] . A theorem, two definitions and a process of modifying a boundary value 
problem that are essential to the applications of this paper, will be recalled in ^1] in order 
to make this paper self contained. The rest of this paper is organized as follows. In $2] we 
prove Theorem 10.61 In §3] we prove a topological index lemma which is a generalization of 
[21 Theorem 1] and [201 Theorem 2]. This lemma may be regarded as a discrete version of 
the Hopf-Poincare Index Theorem with additional terms arising from boundary data. In §U 
we provide several other low complexity examples and their associated propelled surfaces. 
Some of these cases may be analyzed without using the topological lemma of §3j however, 
this lemma provides a uniform approach, hence we will apply it. Finally, $5] is devoted to 
the proof of Theorem 10. 21 

Since one may view a Dirichlet boundary value problem as a special case of a Dirchlet- 
Neumann boundary problem, the main result in [IB] follows from Theorem 10. 21 of this paper. 
However, there is an essential difference between the methods of this paper and those in [16J. 
In some sense the construction of the image surfaces (Sq) in this paper is considerably less 
explicit than the construction in [IB]. At this moment, we are unable to provide a structure 
theory for level curves of the DN-BVP. The structure theorem ([161 Theorem 2.34]) for level 
curves of D-BVP allowed us to cut the domain along the hierarchy of singular level curves 
until we obtain simple regions and then glue back. Thus, in this paper we are led to work 
with the subdomains fij, which turns out to be sufficient for the applications of this paper. 
One may view our results and techniques in this paper as well as in [16] , as providing purely 
combinatorial-toplogical analogues to classical counterparts. See in particular [U Theorem 
4-5], where a connection (in the smooth category) between extremal length of a family of 
curves and Dirichlet Energy of a boundary value problem is exploited. 

Remark 0.9. The assertions of Theorem 10.61 may (in principle) be obtained by employing 
techniques introduced in the famous paper by Brooks, Smith, Stone and Tutte ([8]), in which 
they study square tilings of rectangles. They define a correspondence between square tilings 
of rectangles and planar multigraphs endowed with two poles, a source and a sink. They view 
the multigraph as a network of resistors in which current is flowing. In their correspondence, 
a vertex corresponds to a connected component of the union of the horizontal edges of the 
squares in the tiling; one edge appears between two such vertices for each square whose 
horizontal edges lie in the corresponding connected components. Their approach is based on 
Kirckhhoff 's circuit laws that are widely used in the field of electrical engineering. We found 
the sketch of the proof of Theorem 10.41 given in [8] hard to follow. For a summary of other 
proofs of Theorem I0.6[ a bit of the history of this problem, and generalizations, see Remark 
0.5 in [IB] (as well as [TJ , [T2] , [HI] , and [E]). We include our proof of Theorem 10.61 which is 
guided by similar principles to some of the ones mentioned above, yet significantly different 
in a few points, in order to make this paper self-contained. In addition, the important work 
of Bendito, Carmona and Encinas (see for example [I] , [5j , [6] ) on boundary value problems 
on graphs allows us to use a unified framework to more general problems. Their work is 
essential to our applications and we have used parts of it quite frequently in [16] . and this 
paper as well as its sequel [IT] . 
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I. A REMINDER 

Finite networks. In this paragraph we will mostly be using the notation of Section 2 in [3]. 
Let T = (V, E, c) be a planar finite network, that is a planar, simple, and finite connected 
graph with vertex set V and edge set E, where each edge (x, y) G E is assigned a conductance 
c(x, y) = c(y, x) > 0. Let V(V) denote the set of non-negative functions on V. Given F C V 
we denote by F c its complement in V. Set V(F) = {u G V(V) : S(u) C F}, where 
S(u) = {x G V : u(x) ^ 0}. The set 5F = {x G F c : (x, y) G E for some y G F} is called the 
vertex boundary of F. Let F = F [j 5F and let E — {(x, y) G E : x G F}. Given F C V, let 
r(F) = (F, E, c) be the network such that c is the restriction of c to E. We say that x ~ y 
if (x,y) G E. 

The following are discrete analogues of classical notions in continuous potential theory 

Definition 1.1. ([H Section 3]) Let u G V{F). Then for x G F, the function Au(x) = 
J2y~ x c ( x ' V) ( u ( x ) ~ u (y)) 1S called the Laplacian of u at x, (if x G S(F) the neighbors of x 
are taken only from F) and the number 

E(u) = Y,Au(x)u(x)= c(x,y)(u(x)-u(y)) 2 , (1.2) 

is called the Dirichlet energy of u. A function u G V(F) is called harmonic in F C V if 
Am(i) = 0, for all x G F. 

A fundamental property which we will often use is the maximum-minimum principle, 
asserting that if u is harmonic on V" C V, where V is a connected subset of vertices having 
a connected interior, then u attains its maximum and minimum on the boundary of V (see 
[23 Theorem 1.35]). 

For x G S(F), let {yx, yi, ■ ■ ■ , y m } G F be its neighbors enumerated clockwise. The normal 
derivative (see [TT]) of u at a point x G 5F with respect to a set F is 

£(F)(x)= ^ c(x )2 /)Kx)-«(y)). (1.3) 

The following proposition establishes a discrete version of the first classical Green identity. 
It plays a crucial role in the proofs of the main theorems in [151 HE] and is essential to the 
applications of this paper as well as in its sequel [17] . 

Theorem 1.4. (O Prop. 3.1]) (The first Green identity) Let F dV andu,v G V(F). 
Then we have that 

^2 c(x,y)(u(x)-u(y))(v(x)-v(y)) = ^Au(x)v(x)+ ^ g^( F )( x ) v ( x )- C 1 - 5 ) 

(x,y)eE zG-F xe<5(F) 
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The flux-gradient metric. A metric on a finite network is a function p : V — > [0, oo). In 
particular, the length of a path is given by integrating p along the path (see [9] and [13] 
for a different definition). In [161 Definition 1.9] we defined a "metric" which will be used 
throughout this paper. 

Definition 1.6. Let F C V and let / G V(F). The flux- gradient metric is defined by 

p(x) = if xG 5(F). (1.7) 

This definition allows us to define a notion of length to any subset of the vertex boundary 
of F by declaring: 

Length(<5F) = |^^(F)(:r)|. (1.8) 

xeSF 

In the applications of this paper, we will use the second part of the definition in order to 
define length of connected components of level curves of a boundary value solution. In [T5| 
Definition 3.3], we defined a similar metric (^-gradient metric) proving several length-energy 
inequalities. 

Simple modifications of a boundary value problem. We will often need to modify a 
given cellular decomposition as well as the boundary value problem associated with it. The 
need to do this is twofold. First, assume for example, that L is a fixed, simple, closed level 
curve. Since L R is not (generically) a subset of T^°\ Defintion II .61 may not be directly 
employed to provide a notion of length to L. We therefore add vertices and edges according 
the the following procedure. Such new vertices will be called vertices of type I. 

Let Oi, O2 be the two distinct connected components of L in Q with L being the boundary 
of both (these properties follow by employing the Jordan curve theorem). We will call one of 
them, say 0\, an interior domain if all the vertices which belong to it have g- values that are 
smaller than the g-value of L. The other domain will be called the exterior domain. Note 
that by the maximum principle, one of Oi,0 2 must have all of its vertices with g- values 
smaller than L. 

Let e G and x = e D L. For x G" we now have two new edges (x, v) and (u, x). 
We may assume that v E 0\ and u G 2 . We now define conductance constants c(v, x) and 
c(x, u) by 

c(v, x) = ^ u ]^)-9(u)) ^ = c(v,u)(g(u)-g(v)) ^ 

g{v)-g[x) g(u)-g(x) 
By adding to T all the new vertices and edges, as well as the piecewise arcs of L determined 
by the new vertices, we obtain two cellular decompositions, To 1 of 0\ and 7e> 2 of 2 . Also, 
two conductance functions are now defined on the one-skeleton of these cellular decomposi- 
tions by modifying the conductance function for g according to Equation ( 11. 91) (i.e. changes 
are occurring only on new edges). One then follows the arguments preceding [161 Definition 
2.7] and defines a modification of the given boundary value problem, the solution of which 
is easy to control (using the existence and uniqueness Theorems in [3]). Second, it is easy 
to see that Theorem 11.41 may not be directly applied for a modified cellular decomposition 
and the modified boundary value problem defined on it. Informally, the modified graph of 
the network needs to have its boundary components separated enough, in terms of the com- 
binatorial distance, in order for Theorem 11.41 to be applied. In order to circumscribe such 
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cases, we will add enough new vertices along edges and change the conductance constants 
along new edges in the obvious way, i.e. the original solution will still be harmonic at each 
new vertex and will keep its values at the two vertices along the original edge. Such new 
vertices will called type II. 

2. The case of a quadrilateral 

We now describe the structure of the proof of Theorem 10.61 The proof consists of two 
parts. First, we will show that there is a well-defined mapping from into a set of 
(Euclidean) rectangles embedded in the recatngle Sn- The crux of this part is the fact that 
level curves of g have the same induced length (measured with the flux-gradient metric), 
and a simple application of the maximum principle. Second, we will show that the collection 
of these rectangles forms a tiling of S-% with no gaps. The dimensions of S-r and the first 
Green identity (Theorem 11.41) will allow us to end the proof by employing an energy-area 
computation. 

Keeping the notation of the introduction, we let A, B, C, D be the corners of 1Z ordered 
clockwise with A being the left lower corner and with AB = TZ\ c { t , BC = 7Z t0 p, CD = 7^ r ight 
and DA = 7?.bottom being the boundary arcs decomposition of dlZ. 

Proposition 2.1. For each s G [0,k], the associated g -level curve, l s , is simple and parallel 
to AB, i.e. its endpoints lie on BC and DA, respectively, and it does not intersect ABUCD. 

Proof. Harmonicity of g implies that there exists a path in from B to CD along which 
g is strictly increasing. Since g is extended linearly over each value in [0, k] is attained 
(perhaps at some point on an edge of this path). Let so G [0, k] be any such value. The 
assertion of the proposition is certainly true for s = and for s = k. Therefore assume that 
So G (0, k) and that it is attained at some point which we will denote by Vq. By construction 
v is not an endpoint for l SQ unless Vq G &R, and it is clear that Vq G" AB U CD. Extend l so 
from Vq through triangles and quadrilaterals to a line. It follows by the maximum principle 
that l So is simple and it is not a circle. Also, the intersection of l So with each 2-cell is a line 
segment whose intersection with the boundary of this cell consists of exactly two points, or a 
vertex. Since is finite, l SQ is a closed, connected interval, and by construction may have 
its endpoints only in &R. Let P Vo and Q Vo be its endpoints. To finish the proof, we need to 
show that P vo and Q vo do not belong to the same boundary arc of &R. It is clear that none 
of the endpoints can belong to AB U CD, so suppose (without loss of generality) that they 
belong to BC. Let / = l(P Vo , V\, v 2 , ■ ■ ■ , Q vo ) be the path in BC connecting P vo to Q vo , and 
let V SQ be the polygon formed by l so and the arc /. Attach a copy of it, V so , along /. The 
result is a polygonal disc V so all of its boundary vertices having the same g-value, s . 

Let g be the function which is defined on V SQ by letting g = g on V SQ fl (1Z U &TZ) and 
by letting g(v) = g(v) for every v in the attached copy where v G V SQ is the combinatorial 
symmetric "reflection" of v. By changing the conductance constants (only) along edges in I, 
the fact that g is harmonic in 1Z and since 

fV So )M = 0, (2.2) 

for every v G /. It easily follows that g is harmonic in T> so . However, g has constant boundary 
values, hence g must be a constant (by the maximum-minimum principle). This is absurd. 

□ 
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FIGURE 2.3. A quadrilateral and two parallel level curves. 



Remark 2.4. Similarly, it follows by the harmonicity of g that for each s G [0, k], l s is unique. 



One useful consequence of Proposition 12.11 is that for each s G (0, k) the level curve l s 
separates TZ into two quadrilaterals having disjoint interiors. The first has its left boundary 
equal to AB, its right boundary equal to l s , its top boundary being the part of the top 
boundary of TZ connecting B to the endpoint of l s on it, and its bottom boundary being part 
of the bottom boundary of TZ connecting A to the second endpoint of l s . We will denote this 
quadrilateral by 0\ and its complement in TZ by 02. 

The length of a curve with respect to the flux-gradient metric (see Definition II. 8p . which 
lies on the boundary of two regions, may be computed according to each one of them. The 
length will be said to be well-defined if it does not depend on the region chosen for carrying 
the computation. 



(2.6) 



Proposition 2.5. For every s G [0, k], the length of its associated level curve l s with respect 
to the flux-gradient metric, is well defined and is equal to 

vGCD 

Furthermore, the following equality holds 

dg_ 
dn 



E gwoo 



"> = " E (2-7) 

veCD veAB 

Proof. For s = k the first assertion follows from the definition of the flux-gradient metric 
induced by g. Let s be any other value in [0, k), and let l s be its associated level curve. Let 

v s = / s nr {1) , (2.8) 

and define g at each point of the set V s so that g(y) = g(v), for every v G and 
conductance constants on the added edges so that g is harmonic at each v G V s , and 



|%) = f%) = o, 

on on 



(2.9) 
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where £1 and £2 are the endpoints of l s on BC and AD, respectively (the very last mod- 
ifications are needed only if £1 and £2 are not in T^). Note that the set V s comprises 
(generically) vertices of type I (see the last paragraph in Section [TJ). 

We now apply Green's Theorem (Theorem II. 4 p with u = g and v = 1 over the quadrilateral 
02, which is determined by l s , £iC, CD and to obtain 

= £, e o 2 f n (0,){x) = Et &s SiW(t) + E P ,cn I(0 2 )(P) = (2.10) 



and hence that 

1 E £«wwi - 1 E £woo»i - E |^)w = £ 

te/ 3 pgcd peCD pecz) 

which completes the proof of the first assertion. Note that in order to apply Green's Theorem 
we may need to add vertices of type II and change the conductance constants along added 
edges (see the last paragraph in Section [T]). 

By applying Green's Theorem (Theorem II .4p to 1Z one obtains the second assertion (as is 
the case with the last equality in Equation ( 12. lip , both sides of the equation have the same 
value when computed relative to 7V). In particular, this means that the computation of the 
length of l s with respect to the flux-gradient metric does not depend on which one of the 
two quadrilaterals, 0\ or 2 , it is carried. 

□ 

Given a Euclidean rectangle Q — [0, W] x [0, if] embedded in the Euclidean plane, we 
will endow it with the naturally induced coordinates. Its boundary components [0, W] x 
{0},{0} x [0,H], [0,W] x {H} and {W} x [0,H] will be called bottom, left, top and right, 
respectively. Before providing the proof of Theorem 10. 6[ we need a definition which will 
simplify keeping track of the mapping /. 

Definition 2.12. A marker on a Euclidean rectangle is a horizontal closed interval which 
is the isometric image of [a, b] x {£}, for some t G [0, H] and [a, b] C [0, W\ with a < b. The 
marker's leftmost end-point corresponds to (a,t) and its rightmost end-point to (b,t). 

Proof of Theorem 10.61 Let S-ji be a straight Euclidean rectangle with width W = k and 
height 

H = E ( 2 - 13 ) 



Let C = {Li, . . . , Lk] be the level sets of g corresponding to the vertices in arranged 
in descending (/-values order. We add a vertex at each intersection of an edge with an 
Li, i = 1, . . . , k, (which is not already a vertex in 7~^), and if necessary more vertices on 
edges so that any two successive level curves in £ are at combinatorial distance (at least) 
two. As before, the first group of added vertices is of type I and the second is of type II. 

Starting with x\ = C, we order the vertices {xi = C, . . . ,x p = D} in Li(= CD), as well as 
the vertices on any other level curve, in a monotone decreasing order. Let {y\, 1/2, ■ ■ ■ , yt} be 
the type I neighbors of x\ in the new cellular decomposition oriented counterclockwise (which 
will henceforth be assumed to be the ordering of the neighbors of any vertex). We identify x\ 
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with (k,H) in the coordinates mentioned above, and associate markers {m Xltyi , . . . ,m Xlj y t } 
with x\ in the following way. For s = 1, . . . ,t, the length of the marker m xii y s is equal to 
(the constant) g{x\) — g{y s ) and its rightmost end-point is positioned on the right boundary 
of S_a at height 

H-Y, <x l7 y k )(g(xi) - g{vk)). (2.14) 
fc=i 

For each edge e Uj „ = [u, v] with g{u) > g(v), let Q u ^ v be a Euclidean rectangle with width 
equal to g{u) — g{v) and height equal to c(u,v)(g(u) — g(v)). We will identify a Euclidean 
rectangle and its image under an isometry. For s — 1, . . . , t, we position Q X i,y a i n <SU i n such 
a way that its top boundary edge coincides with m xit y a . By construction and the position of 
the markers, 

Qxi,Vs ^Qx 1 ,y s+1 — m x 1 ,y s+1 - (2-15) 

Assume that we have placed markers and rectangles associated to all the vertices up to 
Xk where k < p; let z 1 be the uppermost neighbor of Xk+i and let Q Xk ,v be the lowermost 
rectangle associated with x^. That is, v is the lowermost vertex which is a neighbor of x^ 
(it may of course happen that v — Z\). We now position the marker m x Zl so that it is 
lined with the bottom boundary edge of Q Xk , v , and its rightmost end-point is on the right 
boundary of S"_4 at height which is given by the obvious modification of Equation f !2.14p . We 
continue placing markers and rectangles corresponding to the rest of the neighbors of Xk+i, 
and terminate these steps when k = p. Note that the right boundary of is completely 
covered by the right boundary edges of the rectangles constructed above, where intersections 
between any two of these edges is either a vertex or empty. 




Figure 2.16. Several rectangles in Sn after the completion of the construction. 



For all 1 < n < fc, assume that all the markers corresponding to vertices in L n _i and their 
associated rectangles have been placed as above in such a way that the following conditions, 
which we call consistent, hold. For [w,v] E with g(w) > g(v) and s E [w,v] a vertex 
of type I, the rightmost end-point of the marker m StV coincides with the leftmost end-point 
of the marker m WjS ] moreover, the union of the rectangles Q WtS ,Qs,v tile Q v>w . Informally, if 
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these conditions are met, this will allow us to "continuously extend" rectangles associated 
with edges that cross level curves along these curves, and therefore will show that edges in 
are mapped in one to one fashion (perhaps in several steps) onto a unique rectangle, 

w,v 

We will now describe how to place the markers and rectangles corresponding to the vertices 
of the level set L n , for n > 2. The rightmost end-point of each marker associated with a vertex 
v G L n and any of its neighbors in L n+ i is placed in Sn on the vertical line corresponding 
to g(v) (the actual height on this level curve is computed by a formula which is an easy 
modification of Equation I2.14p . Observe that v is a vertex in some where belongs 

to L n _i. Choose among all such edges the uppermost (viewed from v G L n ). Let [qo,v] be 
this edge and let m q(hV be its marker. Place m v , w , the marker of v which corresponds to 
an edge [u, iw], with w being the uppermost vertex among the neighbors of v in L n+ i, so 
that its rightmost end-point coincides with the leftmost end-point of m qo ^ v . To conclude the 
construction, continue as above, exhausting all the markers emanating from v, and vertices 
in L n . 

By the maximum principle, our construction, and the fact that all level curves have their 
lengths (with respect to the flux-gradient metric) equal to H, it is clear that the union of 
the rectangles is contained in 5^. 

Proposition 2.17. The placement of rectangles associated to the construction of markers 
as described above is consistent. 

Proof. We prove the assertion by induction on the set of level curves. The assertion is 
obviously true for all rectangles associated with markers emanating from L\, since no such 
marker is a continuation of another one. Let sq be a vertex of type I on L n , n > 2. By 
definition, so is connected to a unique vertex v G L n+ i and to a unique vertex wo G L n -i- 
We first consider the case in which so is the only type I vertex on L n . It is easy to check 
that the following system of equations with the two unknowns c(w , s ) and c(s , i>o) nas a 
non-trivial solution 

c(w , so)(g(w ) -g(so)) = c(s ,v )(g(s ) - g{v )), and (2.18) 
c(w ,s )(g(w ) - g(s )) = c(iu , v )(g(w ) - g(v )). 

The unknowns present the conductance constants to be assigned to [w , s ] and [s ,fo], 
respectively, so that the modified DN-BVP solution function g is harmonic at sq. By the 
construction of the rectangles, this implies that the height of Q So ,v is the same as the height 
of Q(wq,sq), so that they can be glued along the appropriate edges. The second equation 
reflects that the height of the rectangles associated to m WQtSo and m SQtVo is equal to the height 
of the rectangle which one would associate to m WQ)VQ (in the case that there were no type 
I vertices on [u> ,t>o])- In other words, the construction of rectangles is consistent and once 
an edge is split by a type I vertex, the two constructed rectangles my be glued along the 
appropriate edges; thus we obtain the same effect as constructing the rectangle associated to 
the original edge. Note that since g(w ) — g(v ) = (g(w ) —g(s )) + (g(s ) — g(v )), matching 
the widths of the above rectangles is not an issue. 

Assume now that s q is the first vertex of type I in L n which is lower than sq. By def- 
inition, s q is connected to a unique vertex w p in L n _i and to a unique vertex vi G L n+ \. 
Let {si, . . . , s q -i} be the vertices in L n between s and s q , and let {wi, . . . , w p -i} be the 
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vertices in L n _i between wo and w p . Let Qi = Q Wo ,s ,s q ,w p be the quadrilateral enclosed by 
[w , s ] U [w p , s q ] U L n _i U L n , and which contains {w , . . . , w p }, and let Q 2 = Q ao ,vo,8 t ,v, be 
the quadrilateral enclosed by [s , v ] U [s q , vi] UL n U L n+1 , and which contains {s , . . . , s q }. 




Figure 2.19. Viewing Qi U Q 2 . 

In order to prove that the consistent conditions hold for all markers and rectangles created 
in this step, it suffices to prove it at s q ; assuming (without loss of generality) that the first 
marker associated with vertices in L n , that was placed in a consistent way, is m so ^ vo . By the 
construction of the markers (see in particular Equation (12.141) suitably adapted) we need to 
prove that 

EL"l g(Ql)K) + ff (QOlow-leftM + U(Ql)top-leftK) = (2-20) 

EHi 1 g(Q 2 )(*) + S(Q 2 )iow-M t (so) + gCQaWwtfo,), 

where the subscripts "low-left" and "top-left" are posted to emphasize that neighbors in the 
expressions are taken from Qi or Q 2 only. It is easy to check that since g is harmonic at each 
Si, i — 0, ... q (as well as elsewhere), and since s and s q are type I vertices, Equation f)2.20p 
holds. 

Proposition 12.171 

We will now finish the proof by showing that the collection of rectangles constructed 
above tiles Sr leaves no gaps. Without loss of generality, suppose that the collection of the 
rectangles described above does not cover a strip of the form [ti,^] x [ho, hi] in S-r, where 
< t\ < t 2 < k and < Hq < h\ < H. Since g is harmonic there exists at least one path 
whose vertices belong to such that the values of g along this path are strictly decreasing 
from k to 0. In particular, the value t 2 is attained at some interior point of an edge or at 
some vertex of this path. By construction, a gap in a g-level curve (i.e. an arc of a level 
curve which is not covered by the left edges of rectangles) never occurs when t 2 is the g- value 
associated to a vertex in the modified cellular decomposition, T. 

Hence, we may assume that t 2 is attained in the interior of an edge. Let L t2 be the 
corresponding level curve. Recall that L t2 is simple and parallel to CD with its two endpoints 
belonging to AD and BC, respectively (as is the case with all other level curves of g; see 
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Proposition 12. ip . We now follow the construction at the beginning of the proof, and let 
{u\,U2, ■ ■ ■ ,u q } be all the new vertices on L t2 , that is we place a vertex at each intersection 
of an edge in with L t2 . By Proposition 12.61 the length of L t2 (with respect to the 
flux-gradient metric) is equal to H. Moreover, this length is equal to 

E§[(0i)(«*), ( 2 - 21 ) 

i—l 

where 0\ is the interior of the rectangle enclosed by AB, (part of) BC, L t2 and (part of) DA 
(see the end of the proof of Proposition 12. 5p . In particular, in principle we may now place 
in a consistent way, markers and rectangles associated to the collection of edges emanating 
from the vertices {u\, u 2 , . . . , u q } so that L t2 is completely covered by the left edges of these 
rectangles. Since g is extended affinely over edges, every value between ho and h\ is attained 
by g. Repeating this argument shows that all level curves are covered by rectangles. Hence 
the collection of rectangles leaves no gaps in £4. 

Using an area argument, we now finish the proof by showing that there is no overlap 
between any two of the rectangles. Let U be the union of all the constructed rectangles. By 
definition, 

Area(W)= ^ c{x,y){g{x) - g{y)){g{x) - g{y)). (2.22) 

Note that the sum appearing in the right-hand side of Equation (12.221) is computed over 
7W, the induced cellular decomposition. A simple computation (using Equation (12.201) ) 
and the fact that the construction is consistent shows that this sum is equal to the one 
taken over [x,y] G . Hence, the right-hand side of this equation is the energy E(g) of 
g (see Definition II. ip . Therefore, by the first Green identity, applied with u = v = g (see 
Theorem ll.4p . the boundary conditions imposed on g, and the dimensions of S-r, we have 

E(g) = Area(W) = Area(^). (2.23) 

Hence, since the union of the rectangles do not leave gaps, and are all contained in S-%, they 
must tile S^. It is also evident that the mapping / constructed is energy preserving. 

Theorem 10.61 



Remark 2.24. In the forthcoming applications and examples of this paper, we will often 
need to use a slight generalization of Theorem 10.61 First, we will need to allow the domain 
to be sliced. That is, a quadrilateral, two of his adjacent vertices that belong to the left 
boundary (or right boundary), are identified, and possibly a finite number of points on the 
right boundary (or left boundary) are also identified (not necessarily to the same point). See 
the next section and in particular Example 14.131 

One considers a sliced quadrilateral as a quotient of a quadrilateral in the obvious way. 
The construction of Theorem 10.61 goes through with the image being a Euclidean rectangle 
under the appropriate quotient. Note that some of the image rectangles are not going to 
be embedded. However, the embedding fails in a controlled way. For a similar situation in 
the case of a planar pair of pants, and corresponding analysis for higher genus cases in the 
setting of D-BVP, see [I5J Section 4 ]. 
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3. An index lemma 



Let Q be a polyhedral surface with (possible empty) boundary dQ . Let / : Q(°> — > RU{0} 
be a function such that any two adjacent vertices are given different values. Let v G Q^ with 
v G' dQ, and let wi,W2, ■ ■ ■ , Wk be its k neighbors enumerated counterclockwise. Following 
[201 Section 3], consider the number of sign changes in the sequence {f(wi) — f{v), /(ii^) — 
f(v), . . . , f(wk) — f{v), f(wi) — f(v)}, which we will denote by SgCj(w). The index of v G Q 
is defined as 



For the applications of this paper we need to consider the situation in which dQ ^ 0. Let 
v G dQ and let qi,q 2 , ■ ■ ■ ,qi be its neighbors in Q enumerated counterclockwise. Consider 
the number of sign changes in the sequence {f(qi) — fiy), /(q^) ~ f(v), . . . , f(qi) — f(v)}, 
which we will keep denoting by SgcJv). The index of v G dQ is defined as 



Definition 3.3. A vertex whose index is different from zero will be called singular; otherwise 
the vertex is regular. A level set which contains at least one singular vertex will be called 
singular; otherwise the level set will be called regular. 

A nice connection between the combinatorics and the topology is provided by the following 
theorem, which may be considered as a discrete Hopf-Poincare Theorem. 

Theorem 3.4. ([21 Theorem 1], [20~| Theorem 2]) (An index formula) Suppose that Q is 
closed, then we have 



Remark 3.6. Note that due to the topological invariance of x{G) once the equation above is 
proved for a triangulated polyhedron, it holds (keeping the same definitions for Sgc^(-) and 
Ind/(-) as well as the assumption on /) for any cellular decomposition of x{G)- Also, while 
the theorem above is stated and proved for a closed polyhedral surface, it is easy to show 
that it holds in the case of a surface with boundary, where there are no singular vertices on 
the boundary (simply by doubling along the boundary). 

We now prove a generalization of Theorem 13.41 which includes the case of singular vertices 
on the boundary as well as the case in which / admits constant values on some arcs of the 
boundary. Some immediate applications of our generalization will be provided in §4] and £J5] 
providing the control we need on the number of critical points as well as their indices. 

Lemma 3.7. Let Q be a bounded, planar, n-connected domain with dfl as its boundary. 
Suppose that Q U dQ is endowed with a cellular decomposition, denoted by T , in which each 
2- cell is a triangle or a quadrilateral. Suppose that I closed and disjoint arcs are specified 
on the outer boundary of dQ, and that m closed and disjoint arcs are specified on the other 
boundary components of dQ. 

Let f : —> R + be a function which satisfies the following: 

(1) max(/) is attained exactly at each vertex in which lies on any of the I arcs, 




(3.2) 




(3.5) 
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(2) min(/) is attained exactly at each vertex in which lies on any of the m arcs, and 

(3) any two adjacent vertices in T^°\ other than the ones in (1) and (2) have different 
f -values. 

Then we have 

^Ind^) + £ lnd f (v) + = X(O). (3-8) 

ven vedn 

Proof. We first collapse each one of the arcs in dQ on which / attains a maximum or a 
minimum value to a single vertex. The resulting planar domain Q' is bounded and re- 
connected. The cellular decomposition T is changed to a new one T' in the following way. 
Any triangle in with two of its vertices having the same / value is turned into a digon. 
Every quadrilateral with two of its vertices having the same / values is turned into a triangle. 
These changes occur (if at all) only at combinatorial distance which is equal to one from 
the arcs on which / attains a constant value. We now collapse all digons and multi-gons 
connecting two vertices (one of which is on dQ') to a single edge connecting these vertices. In 
particular, we have that x(^) = x(^')> an d T' is comprised of triangles and quadrilaterals. 
Furthermore, / attains its maximum on exactly I vertices in the outer boundary of dQ' and 
its minimum on exactly m vertices in the inner boundary of dQ'. Also, any two adjacent 
vertices in T 1 ^ have different /-values. The indices and the number of the singular interior 
vertices as well as the indices and the number of singular vertices that are on dQ' and are 
not (global) maximum or minimum vertices has not changed. 

We now double Q' along its boundary dQ' and obtain a closed polyhedral surface Q of 
genus x{Q) = 2x(fi')- The index of each interior singular vertex is not changed; however, 
their number is doubled. Let v £ T'^ D dQ' be a singular vertex, and denote by v° this 
vertex in the double of Q'. Note that 

= 2Sgp f {v), (3.9) 

and therefore, 

ind/(?j ) = 2iad/(u). (3.10) 

That is, the index of any boundary singular vertex which is not a maximum or a minimum 
is doubled; however, the number of such vertices is not changed. In the double, we also have 
exactly I vertices on which / attains its maximum, and exactly m vertices on which / attains 
its minimum. It is easy to check that 

lnd f (v°) = 1, (3.11) 

whenever v is a maximum or a minimum vertex. Hence, by applying Theorem 13.41 to Q we 
obtain the following equation. 

2^Ind / (i;) + 2Ind /(^) + {l + m) = X {S) = 2 X (Q). (3.12) 

The assertion of the theorem follows immediately. Let t be the total number of endpoints of 
the arcs on which / is constant. Since = | we will often use an equivalent formulation 
of Equation f l3TT2|) . 

□ 
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4. A FEW LOW COMPLEXITY EXAMPLES AND THEIR SURFACES WITH PROPELLORS 

In this section we will employ the index lemma (Lemma 13. 7p and study a few low complex- 
ity examples. While it is possible to analyze some of the examples in this section without 
applying the index lemma, its usage considerably simplifies the analysis. These examples 
pave the way for the understanding of the general case which will be discussed in the next 
section. 

Example 4.1. (A quadrilateral with two boundary arcs on which / is constant) This ex- 
ample was studied in length in £j2j However, it is worth noting that in this case m — 1 and 
t = 4. Hence, the right hand-side of Equation (13. 12[) must be equal to zero. Since the index 
of an interior singular vertex is smaller or equal to —1, and the index of a singular boundary 
vertex is smaller or equal to — |, it follows that in this case, there are no singular vertices. 
This conclusion is consistent with the assertion of Proposition 12.11 

Example 4.2. (An annulus with one outer Neumann arc) Let A be a planar annulus with 
boundary OA = E±U E 2 , where E\ is the outer boundary. Let ot\ be a closed arc in E\ with 
endpoints Q and P. We solve the DN-BVP as described in the introduction. In particular, 
we have m = 2 and t = 2, and hence that the right hand-side of Equation (13.121) is equal 
to — J- Therefore, the only possibility is that there exists only one singular boundary vertex 
which must belong to Ei\a±. We will denote this vertex by u s , and its associated level curve 
by l s . Since the index of u s is equal to — ~, there must be at least two arcs of l s which pass 
through u s . It follows by the maximum principle that there are exactly two arcs, and that l s 
is simple. Moreover, E 2 U l s comprises the boundary of an annulus which we will denote by 
A[ St E 2 - It follows that A is topologically the union (along l s ) of a topological quadrilateral 
in which two adjacent vertices have been identified and an annulus. Such a quadrilateral 
will henceforth be called a sliced quadrilateral. The following lemma will show that this 
decomposition is geometric. 

Let V = {u s , Vi, v 2 , • • • , Vk] € l s be the set of vertices enumerated counterclockwise. Recall 
that some of these vertices are created due to the intersections of edges in with l s (type 
I), while others may belong to 7~ ( -°\ For any type I vertex, we define the conductance along 
the two new edges it induces according to Equation (12.201) . In particular, if we let g denote 
the solution of the DN-BVP on A which has the same boundary data as g, and the same 
conductance constants on edges which do not have type I vertices, then g and g have the 
same values on and g is a linear extension of g at vertices of type I. 

Lemma 4.3. Let g 2 = g\Ai s e . 2 > ^ e solution of the D-BVP defined on Ai s> e 2 > an d let g\ = 
9\ai s e 2 > solution of the DN-BVP defined on the quadrilateral A \ {Ai St E 2 )° ■ Then the 
length of l s measured with respect to the flux-gradient metric of gi is equal to its length 
measured with respect to the flux-gradient metric of g 2 . 

Proof. Since g is harmonic at each vertex in V which is different from u s , and since the 
Neumann derivative of g at u s is zero, we have that 

E J2 c (y^^9(v l )-g(y)) + ^(A)(u s ) = 0. (4.4) 

We now split the neighbors of each vertex in V other than u s into two groups. For each 
% = 1, . . . , k, let Ai s .E 2 {vi) = {vj, . . . ,vj^} be the neighbors of v t which are contained in 
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(Ai S) e 2 )° and let Qi a> E 2 { v i) — { v l^ +1 ■ ■ ■ > v t } be the rest of its neighbors. Let {it], . . . , u p s } 
be the neighbors of u s in (*4.; Si ,e 2 ) and let {u p+1 , . . . , u g s } be the rest of its neighbors. We 
now rewrite Equation (14.41) in the following form: 

( c(y,v)(g{v)-g{y)) + ^)(#(« ~ 9{y)) ) (4-5) 

+ J2c(u s ,u r s )(g(u s )-g(u:))+ c{u s ,u w s ){g{u s ) - g{u™)) = 0. 

r=l w=p+l 

By the definition of the flux-gradient metric, splitting the sum above into two groups and 
taking absolute values, the assertion of the lemma follows. 

□ 





Figure 4.6. A few level curves in Exam- 
ple 4.2. 



Figure 4.7. The associated surface. 



Example 4.8. (An annulus with one inner Neumann arc) The analysis is similar to the 
one in the previous example. Let A be a planar annulus with dA = E 1 U E 2 , where Ei 
is the outer boundary. Let 0i be a closed arc in E2 with endpoints Q and P. We solve 
the DN-BVP as described in the introduction. In particular, we have m = 2 and t — 2, 
and hence that the right hand-side of Equation (13.121) is equal to — \. Therefore, the only 
possibility is that there exists only one singular boundary vertex which belongs to E2 \ (3\- 
We will denote this vertex by v s , and its associated level curve by l s . Since the index of v s 
is equal to — |, there must be at least two arcs of l s which pass through v s . It follows by the 
maximum principle that there are exactly two arcs, and that l s is simple. Moreover, E\ U L s 
comprises the boundary of an annulus which we will denote by Ai St Ei- It follows that A 
is topologically the union (along l s ) of a sliced quadrilateral in which two adjacent vertices 
have been identified to one, v s , and an annulus. Arguing in a similar way to Lemma 14.31 
shows that this decomposition is geometric. The length of l s measured with respect to the 
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flux-gradient metric of the induced D-BVP on Ai S7 Ei, is the same as measured with respect 
to the flux-gradient metric of the induced DN-BVP in Qa u E2 — A \ (Ai Sj Ei)°- 

Remark 4.9. The surface associated with this example is basically obtained by turning the 
previous one upside down. 

Example 4.10. (An annulus with one outer and one inner Neumann arc) The analysis of 
this case relies on the results and principles set forth in the preceding two examples. Let A 
be a planar annulus with dA — E\ U E 2 , where E\ is the outer boundary. Let ot\ be a closed 
arc in Ei with endpoints Q and P, and let 0i be a closed arc in E 2 with endpoints S and T. 
We solve the DN-BVP as described in the introduction. In particular, we have m = 2 and 
t = 4, and hence that the right hand-side of Equation (I3.12p is equal to —1. By using the 
local structure of an interior singular vertex of index —1, the maximum principle or the fact 
that g has different values on the pairs {P, Q} and {S,T}, one can show that an interior 
singular vertex of index —1 cannot occur. Similarly one rules out the case of a boundary 
singular vertex of index — 1. Hence, the only possible way in which equality may hold in 
Equation ( I3.12p is the case in which there exist two singular boundary vertices, each of index 

— \. With some additional work one shows that E\ \ a.\ contains one of these, say u si , and 
Ei \ P\ contains the other, say u S2 . It follows from the maximum principle that u Sl and u S2 
have different g values. In particular, their associated level curves l Sl and l S2 are disjoint. As 
in the preceding two examples, there are exactly two arcs (of the appropriate level curve) 
meeting at a singular vertex. Hence, A is topologically the union of three pieces. The first 
is a sliced quadrilateral whose boundary consists of E 2 and l S2 which will be denoted by 
Qe 2 ,i S2 - The second piece is an annulus whose boundary consists of l S2 and l Sl which will be 
denoted by A\ i„ ■ The third piece is a sliced quadrilateral whose boundary consists of 1 S1 
and Ei which will be denoted by Qex,i s ■ We have 

Q E2 ,i S2 nA lsi , ls2 =l S2 and Q El ,i ei nA lsi , le2 =i sl . (4.11) 

A simple generalization of Lemma 14.31 shows that the gluing is geometric. That is, for 
% = 1,2, the length of l Si measured with respect to the flux-gradient metric by the induced 
D-BVP on Ai t i B , equals the length measured with respect to the flux-gradient metric by the 
induced DN-BVP on Q.E lt i sl and Qe 2> i S2 i respectively (as before, one needs to add vertices 
of type I and type II, if necessary). 

Remark 4.12. The surface associated with this example is basically the "union" of the sur- 
faces in the previous two. 

Example 4.13. (a planar pair of pants with one outer Nemann arc) Let V be a planar pair 
of pants with dV — E\ U E 2 , where Ei is the outer boundary and E 2 = E\ U E 2 is the inner 
boundary. Let cki be a closed arc in E\ with endpoints Q and P. We solve the DN-BVP as 
described in the introduction. In particular, we have m = 3 and t — 2, and hence that the 
right hand-side of Equation f)3.12p is equal to — |. By applying the maximum principle, the 
case in which there exist three singular vertices each having its index equal to — |, and the 
case of two singular boundary vertices, one of which has index —1 and the other has index 

— |, can be easily ruled out. We are left with the possibility that there exists one interior 
singular vertex, u Sl , whose index equals —1 and one boundary singular vertex, u S2 , whose 
index is equal to — |. Let l Sl be the singular level curve passing through u Sl and let l S2 the 
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singular level curve passing through u S2 . Arguing as in the previous examples, there are 
exactly two arcs of the singular curve passing through u S2 and four arcs passing through u si . 

There are two cases to consider. First assume that u si and u S2 have the same g-value. 
This implies (by using the maximum principle) that they lie on the same singular level curve, 
which we will denote by l s . The second case occurs when u Sl and u S2 have different g- values. 
In the first case the topological decomposition of V is the following. The two annuli Ae^,i s 
and A E 2j s which intersect at u si are attached to the sliced quadrilateral Qi s ,Ei along the 
union of their boundaries, l s . Observe that in this case Qi S! Ei has one singular boundary arc, 
Z s , at u si . 





Figure 4.14. The first case of Exam- 
ple 4.13. 



Figure 4.15. The associated surface. 



In the second case the topological decomposition of V is the following. The two annuli 
Ae\,i s and Ae^i s which intersect at u Sl are attached to the (singular) annulus Ai gi) i a2 along 
their common boundary, Z Sl ; the annulus Ai s t i S2 is attached to the sliced quadrilateral Qi s ^ 



via their common boundary l S2 . 
gluing is geometric. 



It can be shown by a generalization of Lemma 14.31 that the 





Figure 4.16. 
ample 4.13. 



The second case of Ex- 



FlGURE 4.17. The associated surface. 



We finish this section with one more example which illustrates some of the combinatorial 
complexity of higher genus cases. We will not provide a complete analysis of this case and 
leave the completion of the details to the reader. 
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Example 4.18. (a planar pair of pants with two outer Nemann arcs) Let V be a planar 
pair of pants with its boundary V = E\ U E2, where E\ is the outer boundary and where 
E2 — E\ U E\ is the inner boundary. Let <x\ be a closed arc in E\ with endpoints Pi and Q\ 
and let a 2 be another closed arc in E\ with endpoints P2 and Q2, respectively; further assume 
that Pi, Qi, P 2 and Q 2 are ordered counterclockwise. We solve the DN-BVP as described in 
the introduction. In particular, we have m = 3 and t — 4, and hence that the right hand-side 
of Equation (13.121) is equal to —2. There are several cases in which two boundary vertices, 
the index of each is equal to —5, and an interior singular vertex of index which is equal to 
— 1 may occur, and we will now describe two of these. 

Let u si G Q1P2 and u S2 G Q2P1 be the two singular boundary vertices and let tt& be the 
interior singular boundary vertex. First assume that u Sl and u S2 attain the same g values 
and belong to the same level curve which is different from the g value attained at Ub- It 
follows that l s , this singular level curve, is a closed (piecewise linear) curve. It follows by the 
maximum principle that E2 is contained in the domain bounded by l s . Also, If,, the singular 
level curve which passes through w 6 is a piecewise figure eight curve, and (necessarily) the 
g- value of is smaller than that of the g- value of u si . Hence, in this case V has the following 
topological decomposition. A quadrilateral <2 r i g ht which has u si , Qi, Pi and u S2 as its corners. 
A quadrilateral Qi c f t which has u sl ,P 2 ,Q2 and u S2 as its corners. A singular annulus Ai s j b 
with its singular boundary curve being /&, the other boundary curve being l s . The singular 
annulus Ai at i b is attached to Q r i g ht along the right arc of l s , the one which connects u Sl to 
u S2 , and to Q r i g ht along the left arc of l s which connects u Sl to u S2 . 




Figure 4.19. The first case of Exam- Figure 4.20. The associated surface, 
pie 4.13. 
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The singular curve If, = l\ c f t U u& U / r i g ht bounds two annuli, Ae^m* which has as its 
boundary E\ and Ub U Zi e ft, and Asf ,* ri ht which has as its boundary E\ and Ub U Z r ight- These 
two annuli intersect (only) at the vertex Ub- 

We now handle the case in which u si and u S2 have the same g- values but belong to different 
level curves. Let 1 S1 be the singular level curve passing through u Sl , and let l S2 the singular 
level curve passing through u S2 . It is easy to check that lb must intersect QiP 2 in two points 
which we will denote by Si and S 2 , respectively, with Si between u si and Qi and with S 2 
between u Sl and P 2 . Similarly, let T\ be the intersection point of lb with Q 2 Pi which is 
between u S2 and Pi, and let T 2 be the intersection point of 4 with Q 2 Pi which is between u S2 
and Q 2 . Also, l Sl is simple and closed, / Sl C\P 2 Qi = {u Sl }, and the region it bounds contains 
E\. Symmetrically, l S2 is simple and closed, l S2 D Q 2 Pi = {u S2 }, and the region it bounds 
contains E\. 

In this case, the topological decomposition of V is the following. A quadrilateral Qi f t 
which has S 2 ,P 2 ,Q 2 and T 2 as its corners. A quadrilateral Q r j g ht which has Si,Qi,Pi and 
T\ as its corners. A sliced quadrilateral <2 top which has S 2 , Si and u si as its corners; it is 
attached to <2 lcft along the arc of determined by S 2 and Ub, and to Q r i g ht along the arc of 
lb determined by S% and m&. A sliced quadrilateral Qbottom which has Ti,T 2 and u S2 as its 
corners; it is attached to Qi e ft along the arc of lb connecting T 2 and u&, and to Q r ight along 
the arc of lb connecting T\ and Ub- The last two pieces are two annuli, Ae*,i s and A E ^i S2 
that are attached to the above sliced quadrilaterals along 1 S1 and l S2 , respectively. The two 
sliced quadrilaterals intersects (only) at the vertex it&. As before, an extension of Lemma [4731 
shows that in both cases the gluing is geometric. 




Figure 4.21. The second case of Ex- Figure 4.22. The associated surface, 
ample 4.13. 
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5. THE GENERAL CASE - AN m-CONNECTED BOUNDED PLANAR REGION, m > 2 

Proof of Theorem 10.21 Let {0,Pi,P2, ■ ■ ■ ,Pn-i, k} be the set of values of g at the singular 
vertices arranged in an increasing order. We first construct a topological decomposition of 
Q. For i = 0, . . . , k, consider the sub- domain of Q defined by 

fii = {x e Q | pi < g(x) < pi+i}, (5.1) 

(where the value at x which is not a vertex is defined by the affine extension of g). In general 
Qi is multi-connected and (by definition) contains no singular vertices in its interior. Let 
9i — g\ui be the restriction of g to Qi U dQ{. The definition of Qi involves (as in the proof 
of Theorem 10. 6f) introducing new vertices (of type I and type II), and new edges and their 
conductance constants. In particular, each gi is the solution of a D-BVP or a DN-BVP, on 
each one of the components of Q{. 

By applying Equation ( 13 . 1 2 j) in the proof of Lemma 13.71 to g^ and each component of f2; 
whose boundary is a Jordan curve which contains no singular vertices, we obtain that there 
are two cases to consider. First, a component of f2j is simply connected and therefore m — 1 
and t = 4, hence it is a quadrilateral. Second, a component of f2j is not simply-connected. 
In this case we must have m = 2 and t = 0, hence this component must be an annulus. 

We now treat the remaining cases. First assume that the boundary of a component f2^ ordan 
of Qi is a Jordan curve and contains at least one singular vertex which we will denote by 
v s . It follows that the value of g at v s is either pi or Pi+i, and that it does not belong to 
oiU...Ua^U(£ 2 \(/3iU.../3 m )). 

According to the index of v s G Q with respect to g, we now replace a small neighborhood 
of v s in Q by several disjoint piecewise linear wedges. Each wedge has a copy of v s as a single 
vertex, and two consecutive arcs of the associated singular level curve, that are contained 
in the neighborhood, meeting at v s . If v s € dfl, then the only difference from the above, is 
that exactly two of the wedges will contain as one their arcs part of dQ. 



\ Split 




Figure 5.2. Splitting at u 



It follows that after finitely many steps, the boundary of f2 i Jordan is turned into a Jordan 
domain with no singular vertices on it. Hence Lemma 13.71 may be applied to the induced 
Jordan domain and allows us to deduce that f2^ ordan is either an annulus or a quadrilateral. 
It now follows that before the splitting at the singular vertices occurred, f2J° rdan was either a 
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quadrilateral in which two adjacent vertices have been identified or an annulus. One should 
note that singular vertices along the level curves are basically "straightened" along this 
process in such a way that they become non-singular viewed from g^ and £^ ordan with this 
boundary modified to a Jordan curve. 

We must also consider the case in which the boundary of a component fails to be a Jordan 
curve. As Example 14.21 shows, this may already occur in the case of Q being an annulus. 
This case is treated similarly to the previous one we discussed above (see also Figure I5T2]) . 

Thus, we conclude that we may decompose Q into a union (with disjoint interiors) of 
annuli, quadrilaterals or sliced quadrilaterals. We continue the proof by showing that the 
gluing is geometric, i.e. that with respect to the boundary value problems induced on each 
of the components, the common boundary has the same length measure with respect to the 
flux- gradient metric. 

Lemma 5.3. Let L be a connected arc which is contained in TidTj, where Tj is a component 
of Qi and Tj is a component ofQj, for some i and j. Then, the length of L measured with 
respect to the flux-gradient metric induced by g^^ is equal to its length measured with respect 
to the flux-gradient metric induced by gj\rj- 

Proof. The proof is a direct generalization of Lemma 14.31 and follows by applying the 
index lemma (Lemma 13. 7p to rule out several cases. Hence, we will only give the details in 
a few cases. 

Observe that the cases in which L is a Neumann arc, or contains a Neumann arc are 
clearly not possible. By applying the index lemma it can be shown that the cases in which, 
both Tj and Tj are both quadrilaterals is not possible unless L is, without loss of generality, 
the right boundary of Tj as well as the left boundary of Tj (this means that j = % + 1). One 
then uses the fact that the induced D-BVP on Tj U T i+1 is harmonic on L to deduce the 
assertion. One treats the case in which both Tj and Tj are annuli in a similar way; deducing 
that, without loss of generality, j = % + 1, and that the outer boundary component of Tj is 
equal to the boundary component of Tj + i which corresponds to the g value % + 1. Again, 
one uses the harmonicity of the induced D-BVP solution (defined on Tj UTj) on L to obtain 
the assertion. 




Figure 5.4. Viewing Tj and Tj. 
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Assume that (without loss of generality) r, is a quadrilateral and that Tj is an annulus. 
Assume that L is contained in the component of dTj, denoted by dT^ 1+1 , which corresponds 
to the g value Pi+±. Let <9Tj 8 be the second component of Tj (which corresponds to the 

g- value p^), and let L\ +l = <9Tj l+1 \ L. Let P and Q be the endpoints of L. Assume that P 
and Q correspond to the g value Pi + i and are singular vertices. Since L must lie on dTj, it 
follows that the domain bounded by Tj U Tj \ L is an annulus. However, since in this case 
t = 2, we much have (by the index lemma) a boundary singular vertex, this is absurd. 

We now treat one case in which the index lemma does not provide an obstruction for an 
intersection (see §H]for more). The setting is as in the above case, with Tj being this time 
a sliced quadrilateral. Let , A, C, D, B be the vertices of Tj arranged clockwise and let A, B 
be the vertices which are identified. The value of gi on the arc AB is Pi+\. 




Figure 5.5. Viewing Ti and Tj. 

We need to prove that 

By applying Green's theorem to Tj, Tj and Tj U Tj, respectively, we obtain the following 
equations: 

x£BA x£CP x€L x£QD 

" = E tf( T ,>M + E ^raw + Etfraw. < 5 - 8 > 
= E t {T ' ){x) + E f mm + E I*™*) + (5 - 9) 
E §f < r ->M + E ^raw- 

xeQD xedT^ 
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By subtracting the third equation from the first and adding the second equation, Equa- 
tion ( 15. 6 p follows. One verifies by using the method above, that other (finitely many) possible 
cases, lead as well to assertion of the Proposition. 

Lemma 15.31 

We now successively apply the assertions of Theorem 10.61 and Theorem 10.41 to the appro- 
priate components of Q. One needs only observe that the tilling thus obtained is consistent 
as defined in the discussion preceding Proposition 12.171 (see also [IS]). This follows by a 
straightforward generalization of the arguments given in the proof of Proposition 12.171 (see 
the proof of Theorem 10.41 and the proofs of [161 Theorem 0.1, Theorem 0.4]). 

Theorem EH 



Remark 5.10. The analysis of the cone singularities is almost identical to the one carried in 
[I~6| Subsection 4.2]. One observes the following additional cases. The presence of propellors 
results in the creation of new cone singularities of angle tt/2 at each vertex. Hence, under the 
doubling, if such vertex belongs to a unique propellor, the cone angle will change to 7r (see 
for instance vertices Pi, Qi, P2, Q2 in Example 14. 18p . A similar analysis holds if such a vertex 
belongs to two rectangles and a non-singular component of a Euclidean cylinder (yielding a 
cone angle of 4ir in the double). Finally, at the singular vertex of a sliced rectangle the cone 
angle is ir, and the analysis of the changes of this angle under doubling is easy to carry (see 
for instance vertex v s in Example 14.21) . 

Remark 5.11. There is a technical difficulty in our construction if some pair of adjacent 
vertices of has the same g- value (the first occurrence is in Equation (13.1 1) ). One may 
generalize the definitions and the index formula to allow rectangles of area zero, as one 
solution. For a discussion of this approach and others see [191 Section 5]. Experimental 
evidence shows that when the cell decomposition is complicated enough, even when the 
conductance function is identically equal to 1 and the cells are triangles, such equality rarely 
happens (for D-BVP). 

Remark 5.12. The existence of singular curves for g results in the fact that some rectangles are 
not embedded in the target. This is evident by Remark 12.241 and the proof of Theorem 10.21 
Since some of the cylinders or sliced quadrilaterals constructed have a singular boundary 
component, it is clear that some points in different rectangles that lie on this level curve 
will map to the same point. However, this occurs only in the situation described above, and 
since this fact is not of essential interest to us, we will not go into more details. 



26 



SA'AR HERSONSKY 



References 

[1] L. V. Ahlfors, Conformal invariants- Topics in Geometric Function Theory. McGraw-Hill Book Com- 
pany, 1973. 

[2] T. Banchoff, Critical points and curvature for embedded polyhedra, J. Differential Geometry, 1, (1967), 
245-256. 

[3] E. Bendito, A. Carmona, A.M. Encinas, Solving boundary value problems on networks using eguilibrium 
measures, J. of Func. Analysis, 171 (2000), 155-176. 

[4] E. Bendito, A. Carmona, A.M. Encinas, Shortest Paths in Distance-regular Graphs, Europ. J. Combi- 
natorics, 21 (2000), 153-166. 

[5] E. Bendito, A. Carmona, A.M. Encinas, Equilibrium measure, Poisson Kernel and Effective Resistance 
on Networks, Dc Gruyter. Proceeding in Mathematics, (V. Kaimanovich, K. Schmidt, W. Wocss ed.), 
174 (2003), 363-376. 

[6] E. Bendito, A. Carmona, A.M. Encinas, Difference schemes on uniform grids performed by general 

discrete operators, Applied Numerical Mathematics, 50 (2004), 343-370. 
[7] I. Benjamini, O. Schramm, Random walks and harmonic functions on infinite planar graphs using square 

tilings, Ann. Probab. 24 (1996), 1219-1238. 
[8] R.L Brooks, C.A. Smith, A.B. Stone and W.T. Tutte, The dissection of squares into squares, Duke 

Math. J. 7 (1940), 312-340. 
[9] J. W. Cannon, The combinatorial Riemann mapping theorem, Acta Math. 173 (1994), 155-234. 
[10] J.W. Cannon, W.J. Floyd and W.R. Parry, Squaring rectangles: the finite Riemann mapping theorem, 

Contemporary Mathematics, vol. 169, Amer. Math. Soc, Providence, 1994, 133-212. 
[11] F.R. Chung, A. Grigoryan and S.T. Yau, Upper bounds for eigenvalues of the discrete and continuous 

Laplace operators, Adv. Math. 117 (1996), 165-178. 
[12] M. Dehn, Zerlegung ovn Rechtecke in Rechtecken, Mathcmatischc Annalen, 57, (1903), 144-167. 
[13] R. Duffin, The extremal length of a network, J. Math. Anal. Appl. 5 1962, 200-215. 
[14] B. Fuglede, On the theory of potentials in locally compact spaces, Acta. Math. 103 (1960), 139-215. 
[15] S. Hersonsky, Energy and length in a topological planar quadrilateral, Euro. Jour, of Combinatorics 29 
(2008), 208-217. 

[16] S. Hersonsky, Boundary Value Problems on Planar Graphs and Flat Surfaces with Integer Cone singu- 
larities I; The Dirichlet problem, submitted for publication. 

[17] S. Hersonsky, Boundary Value Problems on Planar Graphs and Flat Surfaces with Integer Cone singu- 
larities III; , in preparation. 

[18] P. Hubert, H. Masur, T. Schmidt and A. Zorich, Problems on billiards, flat surfaces and translation 
surfaces, Problems on mapping class groups and related topics, 233-243, Proc. Sympos. Pure Math., 
74, Amer. Math. Soc, Providence, RI, 2006. 

[19] R. Kenyon, Tilings and discrete Dirichlet problems, Israel J. Math. 105 (1998), 61-84. 

[20] F. Lazarus and A. Verroust, Level Set Diagrams of Polyhedral Objects, ACM Symposium on Solid and 
Physical Modeling, Ann Arbor, Michigan, (1999), 130-140. 

[21] H. Masur, Ergodic Theory of Translation surfaces, Handbook of Dynamical Systems, Vol. IB, 527-547, 
Elsevier B. V., Amsterdam, 2006. 

[22] O. Schramm, Square tilings with prescribed combinatorics, Israel Jour, of Math. 84 (1993), 97-118. 

[23] , P.M. Soardi, Potential theory on infinite networks, Lecture Notes in Mathematics, 1590, Springer- 
Vcrlag Berlin Heidelberg 1994. 

[24] A. Zorich, Flat surfaces, Frontiers in number theory, physics, and geometry I, 437-583, Springer, Berlin, 
2006. 

Department of Mathematics, University of Georgia, Athens, GA 30602 
URL: http : //www .math. uga. edu/~saarh 
E-mail address: saarhOmath . uga . edu 



